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Abstract 

We compare the burst distribution of the new (2B) BATSE catalogue to a cosmolog- 
ical distribution. We find that the distribution is insensitive to cosmological parameters 
such as Q and A and to the width of the bursts luminosity function. The maximal 
red shift of the long bursts is ~ 2.1 (assuming no evolution) while Zmax{long) of the 
short bursts is significantly lower Zmax{short) ~ 0.5. In agreement with this relatively 
nearby origin of the short burst we find an indication that these bursts are correlated 
(> 2(7 level at 10°) with Abell clusters. This is the first known correlation of the bursts 
with any other astrophysical population and if confirmed by further observations it will 
provides additional evidence for the cosmological origin of those bursts. 

Subject headings: gamma-rays: bursts — Cosmology 

1 INTRODUCTION 

The recent observation of the BATSE experiment on the COMPTON-GRO observatory 
have demonstrated, quite convincingly, that 7-ray bursts (GRBs) originate from cosmological 
sources (Meegan, 1992). The distribution of the count rate of the bursts in the BATSE IB 
catalogue agrees well a cosmological distribution (Piran, 1992; Mao & Paczyhski, 1992; 
Dermer, 1992; Schmidt, 1992; Wickramasinghe et al., 1993, Mao, Narayan & Piran, 1994). 
Additionally there is a preliminary evidence for the predicted (Piran, 1992, Paczyhski, 1992) 
correlations between the duration, the strength and the hardness of the bursts (Norris, et 
ai, 1994). With the release of the new (2B) BATSE catalogue (Meegan, et a/.1994) we turn 
here to a new analysis of the data, which provides some new evidence for the cosmological 
origin. 

Last year Kouveliotou et a/.(1993) have shown that GRBs are bimodal in duration with 
short (< 2 s) and long (> 2 s) bursts. Lamb, Graziani & Smith (1993) and Mao et a/.(1994) 
have shown that each sub population agrees with a cosmological distribution. Mao, Narayan 
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& Piran (1994) have shown that while BATSE's 1024ms channel is more sensitive to the long 
bursts and the 64ms channel is more sensitive to the short bursts, overall BATSE is more 
sensitive to the long burst which are detected to larger distances. At the same time both 
the long and the short bursts display a surprisingly similar instantaneous peak luminosity 
(to within a factor ~ 2). 

We calculate the likelihood function and show that the observed count distribution in 
the 2B catalogue agrees with a cosmological distribution. The existence of the two sub- 
classes (with different sensitivities) is an important factor that we include in this analysis. 
It would be meaningless to combine both sub-classes to one distribution we treat each sub- 
class separately. Our results support the previous conclusion that the two sub-classes are 
observed to different distances. In the second part of the paper we calculate the cross- 
correlation between GRBs and Abell clusters. We find (at ~ 2 — 3(7 level) a cross-correlation 
for the whole GRBs distribution with Abell clusters and a similar cross-correlation between 
the sub-groups of strong GRBs and of short GRBs (both sub-groups are expected to be 
nearer to us). 

2 Basic Model 

We review here, briefiy, the basic features of a cosmological source count distribution. 
Details can be found in (e.g. Weinberg, 1973) for a general population of sources and in (e.g. 
Piran, 1992) for GRBs. We consider a detector with with a fixed energy range, AD that 
operates for a fixed time span. At, and a source with an energy spectrum: N[E) oc E~". 
The count rate at the detector from a source at a redshift z is: 

where L depends on the luminosity of the source, L, in the relevant energy range, AE", on 
the average energy E and on the observation time, AT: L = L[AE)At/E. Apart from 
normalization which does not vary from one source to another L oc L and hence we will 
delete the" in the following. di[z) is the luminosity distance and for a given 2; it is a function 
of and A. 

Given a population of sources with a luminosity function n[L) and with a rate of events 
n[z) (that is at redshift z we have n[z) events per unit proper volume and unit proper time). 
The number, A^(> C), of observed events with an observer count rate larger than C is: 

[■CO fziCjL) n( z\ 

N{> C) = / niE)dL / ttVt^"' (2) 
JO JO 1 1 -|- 2;) 
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where z[C^L) is obtained by inverting Eq. ??. This theoretical N[C) distribution depends 
on the cosmological parameters: and A, on the source parameters: a, n[z) and n[L). 
Following Schafer et a/.(1992) we assume that all the sources have a spectral index of —1.5. 
This is a crude approximation, but lacking specific spectral indices for the individual bursts 
this is the best that can be done at present. n[z) and n[L) are functions which for simplicity 
we characterize in a few simple ways using one or two parameters. We will characterize 
the distributions in terms of z^ax^ the maximal 2;, from which the detector can detect the 
sources. For a standard candle distribution, n[L) = 5[L — Lq), z^ax is given by: Cmin = 
ivo(l + Zfnax)~"'^'^ /4:'n'df[Zfnax)- Otherwise it is defined as an average over the distribution, 
for example Zjnax{{L)) ■ In feature of the combination of the detector the 

source population and the real source distribution does not end at Zraax- 

3 Maximal Likelihood analysis 

The 2B catalog available in the public domain contains a list of all gamma-ray bursts 
that triggered the BASTF detectors up to March 9th 1993. The catalog contains the ratio 
Cmax/Cmin for 412 bursts, whcrc Cmax IS the maximum count rate and Cmin is the detection 
threshold and Duration table for 434 bursts (Meegan et al. 1994). For a variety of reasons, 
BASTF has a variable background as a function of time, so that the detection threshold 
Cmin does uot remain constant. In order to have a more uniform sample we select a constant 
threshold Ccut (C'cMt( 1024ms) = 287 and Cc„f(64ms) = 72) , and prune the data to contain 
only bursts which satisfy Cmin ^ Ccut and Cmax ^ Ccut (Mao, Narayan & Piran 1993). 
Following Kouveliotou et al. (1993) we define "short bursts" as having Stgo < 2s and "long 
bursts" as having ^tgo > 2s. In order not to exclude bursts which don't have a duration 
record we have used the relation, based on the behavior of 64msec counting and 1024msec 
counting (Mao, Narayan & Piran 1993), for long bursts [Cmax /Ccut) 1024/ [Cmax /Ccut) 64 > 1 
and the opposite relation for the short bursts. Using this relation we add 15 more bursts 
to our set. After excluding bursts which don't satisfy with all the previous requirements we 
have 238 long bursts, and 50 short bursts. We have calculated the likelihood function for 
the multidimensional phase space discussed before. In the following discussion we describe 
several two dimensional cuts through this phase space. 

3.1 Insensitivity to the cosmological parameters: and A 

Figure 1 depicts the contourlines of likelihood function for the long bursts distribution 
in the (il^Zmax) plane. The likelihood function is insensitive to in the range (0.1-1). This 
follows from the relatively modest Zmax value of the currently observed bursts. At these values 
the main effect of on Fq. ?? is to change the length scale. It hardly changes the overall 
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shape of the N[C) distribution. Similarly we find (figure not shown) that the likelihood 
function is insensitive to a cosmological constant A (in the range 0-0.9). Following these 
results we present in the rest of the paper representative likelihood functions for cosmologies 
with = 1 and A = 0. 

3.2 Source Evolution 

The cosmological evolution could be masked by an intrinsic evolution of the sources. 
To examine the effect of intrinsic evolution of the number count we have introduced a z 
dependent local event rate: n{z) = no{l + z)~^ . We describe, in Figure 2, the likelihood 
function in the {zmax^fi) plane. A distribution with no intrinsic z dependence is equivalent 
to a distribution with an intrinsic z dependence with a lower Zraax- This result confirms the 
earlier results of Piran (1992) which were based on the {VjVmax) of the IB catalogue (Note 
that Piran (1992) describes the intrinsic evolution is described as nit) = ngt^' and for low 
Zmax we have (3' —3/3/2). 

3.3 Luminosity Function 

To examine the possible effects of a luminosity function consider a source distribution 
with a luminosity function of the form n[L) = no[L/ Lo)~^ for L > Lq (and N[L) = 
otherwise). Figure 3 depicts the likelihood function in the (s, Zjnax{{L)) plane. The likelihood 
function is quite insensitive to s as long as s 2.2 (note that (L) diverges for s — > 2). This 
shows that, in contrast to previous claims (Mao & Paczyhski, 1992), a narrow luminosity 
function is not needed in order to obtain the observed GRB distribution and GRBs are not 
necessarily standard candles. The insensitivity to the width of the luminosity function was 
confirmed when we tried several other luminosity functions, including a combination of two 
standard candles, a broken power law and a rising power law with an upper cutoff. 

3.4 Long vs short sub-populations 

Figure 2 also depicts the likelihood function for the short bursts sub-population. In 
agreement with Mao, Narayan & Piran (1994) we find that the short burst sub-class is nearer 
than the long burst sub-population. With no source evolution we find that Zmax{-short) < 
1.6 at a 1% confidence level with a maximal likelihood around Zjnax{short) = 0.4. The 
likelihood goes down to 20% of the maximum for 2; 0, hence a nearby (non-cosmological) 
population cannot be ruled out. This is probably due to the relatively low number of bursts 
(50) in this sub-group. For the long bursts we obtain: 1.4 < Zmax{long) < 3.1 at a 1% 
confidence level and a maximal likelihood around Zmax{long) sa 2.1. If we recall that the 
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1024ms channel, that is detecting the long bursts is more sensitive than the 64ms channel 
that is detecting the short bursts we find (again in agreement with Mao, Narayan & Piran 
1994) that both populations have a comparable maximal luminosity. 



4 Cross correlation with Rich Clusters 



We did not find any significant auto-correlation of GRBs in the 2B catalogue. This 
agrees with previous results (Hartman and Blumenthal, 1989) Since GRBs are distributed 
over cosmological distances [z ^ 1) one could not expect a cross correlation with galaxies 
which are relatively nearby. Indeed we did not find any correlation with IRAS galaxies or 
with galaxies in any other optical catalogue. However, Abell clusters (Abell 1958, Abell, 
Corwin & Olwin 1989) are observed up to 2; .15 — .2 and this leads to a possible cross - 
correlation, in particular with strong bursts which presumably originate nearby. We present 
here preliminary results about the cross correlation between the GRBs and rich Abell clus- 
ters. Detailed analysis with respect to various classes of Abell clusters is described elsewhere 
(Kolatt et al., in preparation). 

We carried out the cross correlation with the rich clusters data in the traditional way. 
Let Ngr-d{6) be the number of pairs of GRBs and clusters separated by the angle 6 where for 
the i'th bin 6i_i < 6 < 6i and let Ngr-po{0) be the number of pairs of GRBs and Poissonian 
distributed particles selected according to the same criteria as for the clusters. Then the 
angular cross-correlation function is defined as 

whereas rid and Up^ are the number densities of the clusters and Poisson catalogs respectively. 
In order to account for the different selection criteria we selected GRBs sources with the 
probability 1 — P{S)^ 5 being the terrestrial declination, where 

^ _ # of observing days jS) 

Maximal ^ of observing days 

This procedure prevents us from compensating for the selection (j){S) by multiplying by (j)~^[5) 
and from increasing the noise level in poorly sampled declinations. The selection criteria for 
the clusters were considered by applying the same selection to the Poisson catalog. We 
followed Scaramella et a/.(1991) and used 

MK^) = dex [-A{cosec{\b\) - 1)] X dex [-B{sec{\C\) - 1)] , (5) 

( being the zenithal angle , ( = S — Sfeiescope^ with the appropriate parameters as quoted there 
for the north and south samples. On top of these selection criteria we restricted ourselves 
to Galactic latitude |6| > 30° for both the GRBs and the clusters. 
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A main advantage of the cross correlation method is that it allows to reduce the noise 
level appreciably. In the limit where Upo ^ rid^ the error in the w[6) due to the Poisson 
error is: 

Aw _ 1 

for each 6 bin. The Poissonian errors should still be amplified due to cluster - cluster 
correlations and possible GRBs auto correlation. Figure 4. shows the cross correlation 
between 3616 rich clusters and the GRBs (solid line), the bins' size was determined by the 
GRO angular resolution and was set to 10°. The errors shown for the solid line are not given 
by Eq. 6. Instead, we try here to reject the hypothesis of the GRBs being Poissonianly 
distributed and that the entire cross correlation signal is contributed by the clusters auto 
correlation. We created an ensemble of Poissonian distributed GRBs (same selection and 
number), from which we got the average pair number in each bin as expected for a Poissonian 
model. The error then, was taken to be Aw = \/\ + uj\/< N "pairs >, u being an upper 
limit to the two dimentional J3 integral (u ~ 4). The amplifying factor takes care of the fact 
that the true number of independent pairs is closer to the GRBs - super clusters pairs rather 
than the GRBs - clusters pairs. A more elaborated error analysis is carried out elsewhere 
(Kolatt eA al.^ in preparation). No ensemble average cross correlation between the Poissonian 
model and the clusters was detected. One notices that we can reject the hypothesis of GRBs 
being randomly distributed on the > 2cr level at the 10° bin. 

We proceed by eliminating the presumably non overlapping GRBs and cross correlate 
only the remaining sources. This procedure enables a higher correlation signal level by 
preventing the contamination due to far sources. However, at the same time the noise level 
increases too, due to the reduction in pair numbers. The elimination can be done in two 
ways. We divide the GRBs to strong and weak bursts and compare the cross correlations 
of the two groups with rich Abell clusters. The sum of the two groups does not result to 
the entire sample number since Cmax is not registered for all GRBs. We also divide the 
GRBs to long and short bursts, assuming short bursts are significantly nearer than longer 
bursts. For comparison, we calculate again the cross-correlation function, this time with the 
previously eliminated sources to check consistency of the result with null correlation. Figure 
4 shows the cross correlation of the strong sources (defined by Cmax > Cmaxi'm-edian) for 
the 64ms Channel) with the rich clusters (lower panel). The errors shown, here, are the 
Poissonian errors (Fq. ??) and should be used only as an indication. The actual errors 
derived from a Monte Carlo simulation are larger. The orthogonal set of weak sources is 
shown too. We would like to emphasize that there might be no statistical significant to the 
signal obtained for each bin, but the bifurcation of the two subsets is absolutely consistent 
with the two distinct populations assumption. Figure 4 (upper panel) also shows the same 
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test for sources with short duration (< 2s), again with Poissonian errors (Eq. ??). This 
test turns out again to be consistent with the picture described earher about two clearly 
distinguished populations for the GRBs and with the fact that the short bursts are typically 
nearer than long ones. 

5 Conclusion 

Our results demonstrate that the GRB population is cosmological. Unlike preliminary 
hopes, (Trimble, 1994) we find that the burst distribution is not suitable to serve as a cos- 
mological tool that will distinguish between cosmological models with different cosmological 
parameters. It is possible, but quite unlikely, that much more sensitive detectors that will 
observe the bursts to higher z could do that. This is doubtful since cosmological evolution 
could easily mask the effects of the background geometry. We also find that the distribution 
is rather insensitive to the width of the luminosity function and several luminosity functions 
with different width agree with the data. 

We also confirm the earlier conclusion (Kouveliotou eA ai. 1993) the GRBs distribution 
is divided to long and short bursts. BATSE is more sensitive to the long bursts. If there 
is no burst evolution and the bursts are standard candles then the maximal red shift up to 
which the long bursts are observed is Zmax{long) = 2.lto^7 this corresponds to an event rate of 
S.S^'j'^J- 10~^ events per galaxy per year where we estimated the detection efficiency of BATSE 
as 0.3 and we used a galaxy density of 10~^h^ Mpc""^ (Kirshner et. al, 1983). The rate per 
galaxy is independent of Hq and is only weakly dependent on 0. The typical energy of a burst 
IS 7tl' ■ 10^°F_rergs (where F_7 is the total observed fiuence in 10 ^ ergs/ cm^) for = 1. 
Those numbers vary slightly if the bursts have a wide luminosity function. The rate increases 
and the luminosity decreases if there is a positive evolution of the rate of bursts with time. 
The short bursts are observed only up to a much nearer distances: Zjnax{short) = 0.4"'"^'^ 
(again for standard candles and no source evolution) which corresponds to a comparable rate 
of 6.3~^'^ • 10~^ events per year per galaxy and typical energy of 3"*""^^ • 10'*^F_7ergs. Eor the 
short sub-class, which have only 50 members, there is no significant lower limit on Zraax- 

Recently Norris et a/.(1994) found an increase in the duration of weak bursts by a factor 
of 2 compared to short bursts. They interpret this as evidence to cosmological red shift and 
conclude that the weak bursts originate at 2; 1. Their analysis applies to long bursts 
(longer than 1.5 sec) which, according to our model originate (with standard candles and 
no source evolution) from significantly further distances. This result can be reconciled with 
our analysis if there is a significant source evolution (see Eig. 2) or if there is a sufficiently 
wide luminosity function (the maximal likelihood of the latter case is however only 3% of 
the maximal likelihood of z 2.1). 
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We also found an indication for positive cross correlation of GRBs with rich Abell 
clusters. Abell clusters extend up to 2; .15 — .2 and correspondingly this cross-correlation 
is higher for the stronger half of the bursts (those being nearer) and with the short bursts 
(again those being a nearer sub-population). The cross-correlation is the first association 
of GRBs with any other population of astronomical objects. It supports, of course, the 
hypothesis of cosmological origin of the bursts. It does not mean, however, that the GRBs 
originate necessarily in rich clusters. As it is consistent with the picture in which the bursts 
originate in galaxies which, in turn, are strongly correlated with Abell clusters. This cross- 
correlation indicates that there is some overlap between GRBs and Abell clusters. This is 
self evident for the short bursts for which we found Zjnax{short) 0.4. If the long bursts 
contribute significantly to the cross-correlation this suggests that long GRBs are not observed 
up to Zraax = 2 but to a lowcr Zraax valuc. Just like the results of Norris et a/.(1994) this can 
be reconciled with the N[C) distribution if there is a source evolution which reduces z^ax- 
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FIGURE CAPTIONS 

Figure 1. : Contourlines of hkehhood function in the (0, Zj^ax) plane for standard candles 
sources and with no source evolution. The contours are at likelihood levels 33%, 10%, 
3.3% 1% etc.. of the maximum. The + sign marks the place where the likelihood has a 
maximum. The likelihood is insensitive to 0. 

Figure 2. : Contourlines of likelihood function in the Zj^ax) plane for standard candles 
sources and = 1. (Bottom: Long bursts; Top: short bursts). 

Figure 3. : Contourlines of likelihood function in the [s, Zmax{{L)) plane for = 1 and 
no intrinsic evolution. 

Figure 4. : Cross correlation of GRBs and Abell clusters. Solid line (in both panels) is 
for the entire GRB population, with realistic upper limit errors for a Poissonian model of 
GRBs (see text). The upper panel depicts cross correlation of GRBs for which T50 < 2 ms 
(dotted line). The Lower panel shows the split of the cross correlation when subsets 
according to strength are used. The errors in these two cases are for simple Poissonian 
distribution (Fq. 6) and serve as a lower limit. 
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